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ABSTRACT 

The  atmosphere  causes  errors  in  the  measurement  of  range, 
range  dif ference ,  range  rate,  and  range  difference  rate. 
The  magnitudes  of  these  errors  are  estimated  and  found 
to  exceed  the  maximum  allow  able  errors  for  Mi  strain*  Re¬ 
fractive  index  corrections  must  be  applied  to  the  Mistram 
System  to  reduce  the  errors  to  the  same  magnitude  as  the 
maximum  allowable  errors. 
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1 . 0  Introduction 

The  Mistram  missile  trajectory  measurement  system  derives  its 
Output  from  simultaneous  measurements  of  range  by  a  central  station 
and  remote  stations.  The  output  consists  of  the  missile  range  measured 
at  the  central  station,  Rq  ,  and  range  differences,  DjQ,  determined  from 
the  relative  times  of  arrival  at  the  central  station  and  the  remote  stations. 

The  atmosphere  introduces  errors  in  both  and  DjQ  by  intro¬ 
ducing  changes  in  the  velocity  of  propagation.  The  refraction  errors 
included  in  the  measurement  Rq  range  from  less  than  a  meter  to  a  few 
hundred  meters  depending  upon  the  position  of  the  missile  relative  to 
the  central  station.  The  magnitudes  of  the  range  errors  do  not  appre¬ 
ciably  increase  with  distance  after  the  first  twenty  to  thirty  kilometers 
of  atmosphere  (altitude  above  mean  sea  level)  have  been  traversed. 

This  means  the  relative  accuracy  of  the  range  measurement  increases 
with  range.  The  is  essentially  fixed  by  the  elevation  and  azimuth 
angles  of  the  missile  relative  to  the  site  and  is  approximately  indepen¬ 
dent  of  range.  The  resulting  position  error  caused  by  errors  in  both 
Rq  and  DjQ,  however,  increases  with  increasing  distance. 

The  specified  accuracies  of  the  Mistram  system  are  listed  in 
Table  1.  These  specified  maximum  errors  include  the  average  effects 
of  the  propagation  medium.  These  maximum  error  specifications  hold 
for  elevation  angles,  ao>  between  5°  and  85°  and  ranges  from  55  km  to 
1100  km.  If  the  specified  errors  are  to  be  maintained,  the  errors  due 


to  propagation  through  the  atmosphere  must  be  corrected.  The  expected 
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daily  variations  of  the  average  propagation  errors  also  exceed  the  speci¬ 
fied  maximum  errors.  The  refraction  correction  system  used  with  Mistram 
must  be  able  to  correct  for  the  changes  of  the  atmosphere  with  time. 

TABLE  1 

SPECIFIED  MAXIMUM  ERRORS  FOR  THE  MISTRAM  SYSTEM* 


Range  error 

ARo 

< 

1.  2  x  10"4  km 

Range  difference  error 

ADjo 

< 

9. 2  x  10'6  km 

Range  rate  error 

< 

6.  0  x  10  ^  km/ 

Range  difference  rate  error 

A(  dDjo  1 
dt 

< 

6.  0  x  IQ'7  km/ 

2.0  Range  Measurement  Errors  Due  to  Refraction 

The  Mistram  system  measures  range,  range  difference,  range  rate, 
and  range  difference  rate*  To  make  these  measurements,  the  Mistram 
system  uses  one  central  station  and  from  two  to  four  remote  sites.  The 
range  is  determined  at  the  central  station  by  measuring  the  time  required 
for  a  signal  to  propagate  from  the  central  station  to  the  missile  and  back. 
The  range  differences  are  determined  by  measuring  the  differences  in  the 
times  required  for  the  signal  to  propagate  from  the  missile  to  the  central 
station  and  remote  sites. 

The  Mistram  site  configuration  is  depicted  in  Figure  1.  The  range 
error,  Rq  ,  is  a  function  of  the  position  of  the  target  relative  to  the  central 


*  Page  9,  reference  1. 
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x  -  Mistram  Sites 
o  -  Central  Station 


station  and  the  index  of  refraction  of  the  intervening  atmosphere.  The 
range  error  may  be  approximately  calculated  by  (reference  2): 

10  °  / 

*Ro*!iZZ0J  N(h)dh  in  km.  (1) 

where:  N(h)  =  (n(h)-l)  x  10^,  n  =  index  of  refraction,  and 
a0  =  Elevation  Angle 

The  approximation  holds  for  a©  >10°  and  a  horizontally  stratified  index 
of  refraction.  On  the  average  at  Cape  Canaveral  the  integral 


Io(»)  =  10-6 


N(h)dh  =  2.  5  x  iO'^km. 


and  represents  the  range  error  for  a  vertical  ray  path  through  the  atmos 


phere.  The  average  expected  range  error  is 
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~  htl 

“  sincfe 


2.  5  x  10 
sin  Qo 


km 


2.  5 
sina0 


meters 


(2) 


where  the  height  of  the  target  is  greater  than  30  km. 

For  target  altitudes  less  than  30  km,  the  range  error  will  be  less. 

The  form  of  the  radio  refractivity  profile,  N(h),  is  roughly  ^.n 
exponential  decrease  with  height.  For  a  given  elevation  angle,  the  error 
per  unit  path  length  decreases  in  an  exponential  fashion  with  height.  The 
relative  accuracies  of  the  range  errors  then  increase  with  range  or 
height  for  a  given  elevation  angle.  For  elevation  angles  greater  than 
10°,  the  angle  clq  is  approximately  the  elevation  angle  for  a  straight  line 
between  the  target  and  the  central  station.  For  increasing  ranges  and 
constant  elevation  angles,  the  increase  in  total  range  error  is  negligible 
for  altitudes  greater  than  30  km. 

For  elevation  angles  less  than  10°,  the  range  error  cannot  be  ac¬ 
curately  computed  using  equation  (1).  The  exact  equation  for  range  error, 


assuming  horizontal  stratification,  is  (reference  3) 

rhfc  N(h)  dh  tyt  dh 

^  =  10  /  sin('^.) 

hr 


f  = 


R. 


in(a0  +<£“T')  SR 


where  T  is  the  angular  bending  of  the  ray  path,  0 
is  the  central  angle  (relative  to  the  center  of  the 
earth)  subtended  by  the  central  station  and  the 


target,  and  R  is  the  straight  line  distance  from 
SR 


(3) 


the  central  station  to  the  target* 
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The  ranges  of  values  of  AR$  that  may  be  expected  may  be  estimated  using 
the  NBS  exponential  atmosphere  (reference  4). 

TABLE  2 

RANGE  ERRORS  CALCULATED  USING  THE 
CRPL  EXPONENTIAL  ATMOSPHERE  FOR  TARGET  HEIGHT  OF  425  KM* 


Elev.  Angle 

ARo  for  Ns=200 

ARC  for  Ng^450 

0° 

65  m. 

192  m. 

3°  43 1 

21  m. 

29  m. 

11°  271 

9  m. 

■  10  m. 

A  sampling  of  84  radio  refractivity  profiles  measured  at  Cape  Canaveral 
were  studied  in  detail  by  Thayer  and  Bean  at  NBS.  The  ranges  of  values 
calculated  from  these  profiles  are  presented  in  Table  3. 

The  range  errors  given  in  Tables  2  and  3  are  for  an  atmosphere 
whose  index  of  refraction  is  horizontally  stratified.  Any  real  atmosphere 
is  not  horizontally  stratified  but  includes  tilted  layers  and  local  inhomogen¬ 
eities  such  as  clouds.  The  effects  of  the  local  inhomogeneities  in  the  at¬ 
mosphere  may  be  investigated  by  using  the  stratified  atmosphere  calcula¬ 
tions  and  adding  correction  terms.  The  effects  of  tilted  layers  may  be 
investigated  by  moving  the  origin  of  the  coordinate  system.  The  range 
error  may  be  expressed  by  (Reference  3) 


~  f*  * /'  JVWf;  +  “  (“■  mkdR 

J  sin  a  J  sin  a  k=l  J 

rki 


*  Extracted  from  tables  in  reference  4. 
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~  I  NCT(h)dh 

[T  8iA a'  +  £  (%  -  rk1>  (4) 

th 

where  the  AN^.  is  the  k  deviation  from  the  stratified  atmos¬ 
phere,  Na(h)  refers  to  the  assumed ^stratified' atmosphe re,  and 
-  r^  is  the  path  length  through  the  k**1  local  dis¬ 
turbance  * . 

The  magnitude  of  each  disturbance  is  assumed  small  so  it  does  not  affect 
the  total  amount  of  bending  along  the  path.  For  a  cloud,  the  AN  value  is 
about  10  (reference  5).  If  a  length  of  about  1  km  is  assumed,  an  increase 
in  range  error  of  about  0.01  meters  may  be  expected.  Since  this  is  neg¬ 
ligible.  compared  to  the  total  error  accumulated,  the  effects  of  local  in¬ 
homogeneities  will  only  be  noted  (reference  6). 

TABLE  3 


RANGE  ERROR  CALCULATED  FROM  MEASURED  PROFILES* 


Elev.  Angle 

Ht. 

^min: 

AR 

max, 

0° 

1  km. 

23  m. 

41  m. 

59  m. 

5  km. 

46  m. 

78  m. 

107  m. 

70  km. 

65  m. 

107  m 

145  m. 

5 1 ' 

1  km. 

3.  7  m. 

5.  8  m- 

7.  1  m. 

5  km. 

14.  5  m. 

21.  2  m. 

24.9  m. 

70  km. 

27.  5  m. 

38*  7  m. 

43.  6  m. 

1 1°  27' 

1  km. 

0.  96  m. 

1*49  m. 

1. 84  m. 

5  km. 

4*  05  m. 

5.81  m. 

6.  75  m. 

70  km* 

8.  6  m. 

11.9  m. 

13.  1  m. 

♦Extracted  from  Table  1,  pg.  4  of  Reference  2. 
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The  phenomenon  of  ducting  will  also  increase  range  and  range  dif¬ 
ference  errors.  The  approximate  amount  of  increase  depends  upon  the 
target,  site  geometry  and  the  atmospheric  conditions.  In  general,  ducting 
will  not  be  a  problem  in  the  area  specified  for  the  operation  of  the  Mistram 
system  (references  6  and  7). 

The  refractive  range  errors  vary  between  a  part  of  a  meter  and 
two  hundred  meters  depending  upon  the  atmospheric  conditions  and  the 
target,  site  geometry.  The  variation  in  range  error  may  be  as  large  as 
fifty  percent  of  the  average  range  error  (Table  3).  The  maximum  range 
error  specified  for  the  Mistram  system  is  0,  12  meters.  The  expected 
errors  and  changes  in  errors  due  to  atmospheric  refraction  are  many 
times  the  maximum  allowable  error. 

3.  0  Range  Difference  Measurement  Errors  Due  to  Refraction 

The  refraction  error  included  in  the  measurement  of  Djq  may  be 
computed  if  the  range  errors  for  both  the  central  station  and  the  remote 
station  are  known, 

ADjo  =  (Rj  +  ARj-R0-AR0)  -  (Rj-R0)  =  ARj  -  AR0 

For  elevation  angles  greater  than  10°, 

AD:  =  I  (“)  I — r- -  -  — ^ - 1  (5) 

J°  °  [  sin  aj  sin  aQ  J 

where  the  index  of  refraction  profile  is  assumed  to  be 

horizontally  stratified  and  identical  for  both  ray  paths. 

For  purposes  of  analysis,  the  target,  site  geometry  may  be  reduced 
to  the  plane  figure  depicted  in  Figure  2.  An  effective  baseline,  Bjoe, 
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is  introduced  to  compensate  for  the  azimuth  angle,  shown  in  Figure  1. 
The  maximum  error,  ADjomax'  occurs  when  the  target  is  along  the  ex¬ 
tended  baseline  or  when  =  B;n. 


For  large  Rq  the  ray  paths  are  approximately  parallel*  At  the  target 
the  elevation  of  the  ray  relative  to  the  horizontal  is 

Oto'  =  ao  +  0  -  To 

atraj  +  *+^je-Tj 

B . 

where  0.  =  Joe  ,  r  =  radius  of  the  earth  at  MSL. 


Je 


-  a.  =  0j„  +  +  (ato-atj)  = 


O 

J  ')« 

Using  Equation  (5) 

~  W  Yje  CQt  ao  C8C  ^ 


je 


AD. 


J° 


l  -  Yje  cot  aQ 


For  target  altitudes  greater  than  30  km, 


(7) 


-  T  i  <  A 

o  J  Te 

cl  -  a.-  <  0. 
to  “tj  *je 


As  a  crude  approximation, 


do  *  0L\ 


J  Tje 


=  <k 


je 
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~  Io(“)  0jeC8COo  c°ta0 

ad. = - — : - *  a° >  10 

J°  1  -  0je  cotao 


(8) 


For  elevation  angles  less  than  10u,  the  range  error  may  be  related  to  I^h) 

by  (reference  3) 

AR.  =  — i -  IJh) 

J  sin  aej  ° 

where  the  effective  angle  q^j  is  a  function  of  the  target, 

site  geometry  and  the  radio  refractivity  profile  N(h). 

This  equation  may  be  used  to  define  the  angle  a  .  given  AR*  and  I0(h). 

ej  j 

A 

sin  a  .  =  - 

eJ  ARj 

The  effective  angle  for  one  path  may  be  related  to  the  effective  angle 
for  the  other  path  by 

r  , 

°eo<10U 


°eo  -  aej  *  loW 


Li.  .  jJ  • 

[t*o  ARj  J  ’ 


Since  for  large  R0,  the  effective  angle,  otej  would  bear  the  same  rela¬ 
tionship  to  the  angle  parameters  of  ray  path  R^  as  aeo  does  to  ray  path 
R0,  the  effective  angles  may  be  related  by  the  crude  approximation; 

aeo  '  ^je 


The  range  difference  error  then  is  approximated  by 
^  Io(0D)  foe  CQt  P  CSCP 
1  -  0je  cotP 


AD. 

jo 


(9) 


where  p 


_  fa0,  Oo  >  10? 

l?eo»  ^o  ^  10° 


and  ht  >  30  km 


In  terms  of  I  (®)  and  ARQ,  the  range  difference  error  may  be  expressed 


by 


AD- 

jo 


I°(")  0ie  cot  P 
Sin  3  1  -  0je  COt  |3 
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ARp  0je 
tan  p  -  0jfe 


^Dj°max 


ARp  0j 
sin  g  ^ 


for  g  <  30°  or  530  mr 


ARo0j 
Xo  H 
ARo 


* 


=  (AR0)8 


B 


22- 


rolol00 > 


10. 


(10) 


(11) 


which  holds  for  IQ  (®)  >  >  ARQ  0j  or  in  most  cases. 

AD.q  =*  1.92  ( ARq)2  x  10 -3  meters  for  ARQ  in  meters  ( 1 2) 
^°max  J 

g  <  30°,  a  baseline  of  30 . 48  km  (100  k  ft)  and  I©  (“)  -  2.  5  x  10"  ^km 
(5  year  mean  fox'  Cape  Canaveral), 

The  last  expression  is  a  good  approximation  for  angles  below  30°  . 

For  angles  between  30°  and  85°  it  will  hold  within  an  order  of  magnitude 

(RQ  large)  since,at  85°,  tan  g  is  an  order  of  magnitude  larger  than  sin  g. 

This  expression  therefore  may  be  used  to  determine  the  upper  bound  on 

the  range  difference  errors  expected  for  the  entire  volume  of  space  of 

interest  to  Mistram.  The  maximum  values  being  obtained  for  targets 

along  the  extended  baseline,  Bja 

For  an  elevation  angle  of  200  mr  (11.5°)  the^average  of  I0(»)  for 

Cape  Canaveral  is,  using  (2), 

AR  =■  — - —  ■  -  *  =,  12.  6  meters. 

°  sin(0.  2  rad) 

The  values  of  range  error  tabulated  above  (tables  2,  3)  shows  for  this 


elevation  angle, 
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8.  6m  <ARQ  <  13.  1  meters. 

For  this  case,  using  (9)  and  a  baseline  BjQof  30.  48  km.  ( 100  k.  ft. ) 

ADj0  -  0.  30  meters 

or  using  ( 1 2) 

AD  =  0.  30  meters, 
jo 

The  errors  calculated  are  for  an  elevation  angle  of  11.5°  and  applies  for 
ranges  greater  than  150  km.  At  this  elevation  angle,  the  range  error  varies 
from  about  8  meter 9  to  about  13  meters.  The  corresponding  range  difference 
error  varies  from  0.  12  meters  to  0.  33  meters.  In  both  cases  both  the  aver¬ 
age  expected  errors  and  the  day  to  day  changes  in  the  average  error  far  ex¬ 
ceed  the  maximum  allowed  errors  of  ARQ  =  0.  12  meters  and 

max 

AD  =  0.  009  meters  (table  1). 

J°max 


The  effect  of  local  inhomogeneities  in  the  atmospheric  index  of  refraction 
on  ARq  was  noted  above.  The  same  holds  for  ADjq,  An  inhomogeneity  such 
as  a  cloud  that  lies  along  the  ray  path  between  either  the  remote  site  and 


the  missile  or  between  the  central  station  and  the  missile  and  not  along 


the  other  ray  path  will  cause  a  change  in  ADjq. 

A  change  in  ADjQ  of  about  0.  01  meters  would  result  if  an  inhomogeneity 
occurs  of  the  same  size  and  average  index  change  as  computed  above. 


This  change  in  AD^o  Is  also  of  the  same  order  of  magnitude  as  the  maximum 
allowed  range  difference  errors.  Larger  disturbances  such  as  frontal 


movements  will  cause  even  a  more  drastic  change  in  AD- 


TM-3422 


12. 


4.  0  Range  Rate  and  Range  Difference  Rate  Errors 

Errors  in  range  rate  and  range  difference  rate  will  also  be  caused  by 
atmospheric  refraction.  The  errors  are  caused  by  the  rate  of  change  of  the 
range  error  caused  by  the  atmosphere. 

dRa  d  .  dRSR  d 

if  -  3?  <rsr  +  ®.I  -  ~dt  +  7,  <“«> 

The  range  rate  error  then  is 

A  (  *“)  =  It  (A»o) 

and  is  related  to  both  the  rate  of  change  of  the  target  along  a  particular  ray 
path  (g  =  constant)  and  the  rate  of  change  of  g.  For  motion  along  the  ray 
path  and  a  target  height  greater  than  30  km,  the  range  rate  error  is  negligible. 
The  error  in  this  case  is  negligible  because  the  change  in  range  error  per 
unit  change  in  range  is  negligible  for  hj.  >  30  km.  The  range  rate  error  then 
may  be  considered  due  to  change?  in  g, 

=  -  cote  If  cetera/ sec  (13) 

for  ARq  in  meters  and  ht  >  30  km. 

Similarly,  for  range  difference  rate  errors  due  to  refraction, 

(ADj0)  "=  ^  (1.92  (bRg)3  X  10"3  )  =  -  2AD.q  cot  3  If  meters/ sec 

(14) 

for  ADjQ  in  meters  and  h^  >  30  km. 

The  rate  errors  described  above  kre  for  a  spherically  symmetric 

atmosphere.  These  estimates  do  not  include  the  rate  errors  caused  by 

the  motion  of  the  ray  pkth  through  inhomogeneities  in  the  atmosphere. 


/ 
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5 .  0  The  Effect  of  Range  and  Range  Difference  Errors 

The  range  errors  also  contribute  to  errors  in  target  position.  A 

detailed  analysis  of  position  errors  is  impossible  without  a  lengthy 

consideration  of  the  use  of  the  Mistram  outputs.  The  general  effect  of 

errors  in  Rq  and  Djo  may  be  inferred  by  considering  the  magnitude  of  the 

position  error,  L,  depicted  in  figure  3. 


Figure  3 

The  measured  distances  from  the  central  station  to  the  target, 
xm  and  Ym  are 


2RoDjo+(Djo)2-(Bjo)2 


2  B 


(15) 


y  =  V*  2  *  *  2  ( 

7  m  V  o  m 

where  RQ  is  large,  the  rays  are  assumed  to  propagate 
in  straight  lines,  and  the  target  is  assumed  to  be  above  the  extended  base 
line.  The  position  errors  are  related  to  A  RQ  and  A  Dj0  by 


I  Ax  |  = 

I  Ay  |  = 


AR0  Djo  +  ADj0(RQ  +  Dj0) 


Rq  A  Rq  -  xm  A  x 

ym 


(17) 

(18) 


/ 
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14. 


L  =/(A  x)2  +  (A  y)2  .  (19 

The  magnitudes  of  the  position  error,  L,  for  the  case  computed  above  is 

L  -  15.13  mete r s 

where  |  A  x  |  =  1 4.T  Zm,  |A  y  |  =  5 .  44m 

xm  =  146.  41  km,  ym  =  32.  60  km 


ARq  =  12.  6m, 


A  Djo 

=  0-  3m 

o’ 

II 

30. 48  km 

,  B. 

=  30. 48  km. 

The  asymptotic  behavior  of  the  position  errors  for  large  RQ  may  be 
deduced  from  the  above  expression.  For  large  RQ,  ARQ,  and  Dj0  are  essentially 
constant. 


xm  =  ^  Ro  Dj0  4  (Djo^  ~ 
2  B~ 


Ro  Djo  =  k,  R 


Ym  “  A/Ro  ”  xm^ 


Ro  V1  -(  Djo 


k2  R0 


Ax  |  —  A  R0  Dj0  +  A  Djo  ( Ro  ^  Djo) 


|  A  y|  =  R0  ARo-  xm  A  x 

ym 

L  =  V(Ax)2  +  (A  y)2 


_  ~  R°(4-$)  =  k3R»  (22> 

R  (A  R  -  k,  k,  Rn) 
o  1  o  1  3  O'  Ir,  1 _ 


k2  Ro 

Ro  k3~\/1  YM  2  =  k4  R0 


kl  k3  R0  (23) 
k2  . 

(24) 
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15. 


The  magnitude  of  the  position  error  therefore  asymptotically  approaches  a 

linear  increase  with  RQ.  For  the  case  computed  above,  L  increases  linearly 

with  R  for  R^  >  15,000  km 
o  o 

The  errors  in  range  difference  and  position  have  been  considered  in 
detail  for  the  target  along  the  extended  base  line,  0j6*  =  0>j.  In  general,  the 
target  will  be  at  some  angle  with  respect  to  the  base  line.  The  effective  great 
circle  arc  between  the  remote  site  and  the  control  station  will  be  smaller  or 


From  (10),  the  range  difference  error  approximately  is 
A  Djo  =  (AR0)a  0je' 

1^*)  “  V 


Similarly,  from  (22) 


Therefore, 

L  =  "V(  Ax)a  +  (Ay)2  «  0..- 

‘  'J" 


In  the  extreme  case,  when  the  range  difference  error  is  negligible, 
0je  =  Of  the  magnitude  of  the  position  error  is  due  to  ARq  albrie  and  inde-r 
pendent  of  range  as  ■■5  ,f’  *'*■  'v-'  ’•  ' 
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The  Mistram  system,  however,  requires  data  from  a  remote  site, 
central  station  pair  along  the  other  baseline,  to  estimate  position  in 
three  dimensions.  The  error  from  that  station  will,  in  general,  be  a 


maximum,  0.  *=  0.-  ,  and  the  composite,  three  dimensional  error  will 

je  J 

asymptotically  approach  a  linear  increase  with  R^. 


6.  0  Conclusic 


The  relative  range  errors  and  range  difference  errors  decrease 
with  increasing  distance  from  the  central  station.  The  magnitude  of  the 
position  error,  however,  increases  with  distance  from  the  central  sta¬ 
tion.  The  average  refractive  errors  exceed  the  maximum  allowed  errors 
(Table  1)  by  one  to  three  orders  of  magnitude.  Refractive  index  corrections 
must  be  applied  to  the  Mistram  outputs  to  maintain  refractive  errors  less 
than  the  maximum  allowable  overall  system  errors. 

, (  The  daily  variations  in  the  range  and  range  difference  errors  due 
to  changing  atmospheric  conditions  also  may  exceed  the  maximum  allow¬ 
able  errors  by  one  to  two  orders  of. magnitude .  Errors  due  to  local  in¬ 
homogeneities  may  be  of  the  same  orders  of  magnitude  as  the  specified 
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maximum  allowable  errors.  To  maintain  the  specified  accuracy,  the 
Mistram  system  must  periodically  sample  the  effects  of  the  atmos¬ 
phere  and  correct  the  output  data.  The  refraction  corrections  must 
also  be  applied  to  each  ray  path  separately. 


£ ZrL±AA  & == £ 

R.  K.  Crane 

Electronic  Warfare  Department 
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